Extracting a. from — ► K°K° Decays 
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We propose a new method for obtaining the CP phase a, based on measurements of B% — > K°K°, 
along with theoretical input. Due to the similarities of QCD factorization (QCDf) and perturbative 
QCD (pQCD), this input is basically the same for each of these models. Although the theoretical 
error is large at present, many of the contributing quantities will be better known in the future, 
leading to a smaller error. One outstanding question is the extent to which the precision on quark 
masses, especially m c /mb, can be improved. If one assumes that new physics is not present, the 
method can be used to predict the values of CP asymmetries in B® — » K°K°: 0.02 < A dir + < 
0.125. A result outside this range would signal the existence of long-distance effects beyond those 
included in models of nonleptonic decays based on factorization (such as QCDf and pQCD), or the 
presence of new physics. 
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A great deal of work has gone into finding a large va- 
riety of ways of measuring a, (3 and 7, the three interior 
CP-violating angles of the unitarity triangle [l[ . The idea 
is to find a discrepancy between different ways of mea- 
suring the same CP phase, or with the predictions of the 
standard model (SM). This would indicate the presence 
of physics beyond the SM; it would then be necessary 
to include extra (measurable) new-physics parameters in 
the amplitude of any process to be analyzed Q. 

In the present paper, we present a new method of mea- 
suring a. It involves the decay B^ — ► K°K°, along with 
theoretical input. At present, the theoretical error is 
large. However, it is a function of quantities such as the 
renormalization scale, light-cone distributions, and form 
factors, all of which will become better known with time. 
It is also a function of quark masses. If these can be de- 
termined more precisely, especially m c /mb, this method 
can eventually provide a rather accurate way of extract- 
ing a, and may help to reveal new physics. 

There are several amplitudes contributing to the de- 
cay B® — > K°K°, which at the quark level is 6 — > dss. In 
diagrammatic language 0] , the largest is the penguin di- 
agram pen. This amplitude receives contributions from 
each of the internal quarks u, c and t. However, using 
the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, one can eliminate the c-quark contribution and 
write 



pen = v: b v ud (p u - p c ) + v; b v td (p t - p c 



(i) 



Note that P u and P c arise mainly due to rescattering of 
the tree diagrams. The full amplitude for B® — > K°K° 
can therefore be written 



A = A(B% -» K°K ^ 



v: b v ud [(p u ~ p c ) + 

+ V t * b V td [(P t -P c 



In the Wolfenstein approximation [4| , only the CKM ma- 
trix elements V u b and Vt d are complex. V* b and Vtd con- 
tain the weak phases 7 and — (3, respectively, so that the 
relative weak phase is (3 + 7 = it — a. The penguin 
amplitudes P t contain strong phases. The amplitude A 
describing the conjugate decay B d — > K°K° can be ob- 
tained from the above by changing the signs of the weak 
phases. 

The amplitudes A and A thus depend on four unknown 
parameters: the magnitudes Vtc = \ [{Pt — P c ) + ■ ■ -]V t * b V td \ 
and Vuc = \[(Pu - P c ) + —]V* b Vud\, the relative strong 
phase A = 5 UC — S tc , and a. However, there are only three 
measurements which can be made of B ( j — > K°K°: the 
branching ratio, and the direct and mixing-induced CP- 
violating asymmetries. These yield the three observables 



X = X - (| A\ 2 + \A\ 2 ) =Vl c + Vl - 2V uc Vtc cos A cos a , 

Y = X - (| A\ 2 - \A\ 2 ) = -2V uc V tc sin A sin a , (3) 

Z T = Im (e~ 2ll3 A*A) =V 2 UC sin 2a - 2V uc V tc cos A sin a . 

In the above, 2(3 is the phase of B^-B d mixing. 
It is useful to define a fourth observable: 

Z R = Re {e~ 2lP A*A) 

= V 2 UC cos 2a + V 2 C - 2V uc V t c cos A cos a 



(4) 



The quantity Z R is not independent of the other three 
observables: 



Zl 



X — Y 



Z 



(5) 



■(2) 



Thus, one can obtain Z R from measurements of X, Y 
and Zj, up to a sign ambiguity. 

Note that the three independent observables also de- 
pend on four theoretical parameters (V uc , Pta A, a), so 
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that one cannot obtain CP phase information from these 
measurements [5|. However, one can partially solve the 
equations to obtain 



2 _ .ZhCOS 2a 

' tc 



Zj sin 2a - X 



cos 2a — I 



(0) 



In Ref. [6|], two of us (AD,DL) noted that another con- 
straint on Vtc can be obtained from B s -> K°K°. The 
argument goes as follows. To a good approximation, the 
amplitude for B s — ► K°K° is given by 



K U K 



+ v; b v ts (Pi~Pc) , 



(7) 



where the prime indicates a b — > s transition. The first 
term is much smaller than the second since V* b V us <C 
Vf^Vis, and can be neglected. The decay B s — > AT ^ 
thus basically involves only P' tc = — P' L )V* b V ts \, and 
this quantity can be obtained from its branching ratio. In 
the limit of perfect flavour SU(3) symmetry, P[ c — Vtc 
apart from known CKM matrix elements. Thus, the mea- 
surement of P' tc gives the necessary additional condition 
for extracting a from B d — ► K°K° (a connection with 
Bd — > 7T7T can be found in Q). However, SU(3)-breaking 
effects will affect the relationship between P' tc and Vtc- 
These are 0(30%) and are unknown, leading to a large 
theoretical error. In Ref. [6j this error was reduced to 
0(5%) by considering two decays, e.g. B d — » K a K° and 
B d — > if* A"* , and using double ratios. 

In the present paper, we point out that there is an 
alternative additional piece of information. In order to 
see this, it is simpler to eliminate the i-quark contribution 
from pen 1 writing 



A(B» d -> K U K») = V* b V ud T + V c * b V cd P 



(8) 



where P = (P c - P t ) + ... and T = (P u -P t ) + ... are com- 
plex quantities. (The ... represent smaller contributions.) 
Compared to the previous parametrization [Eq. @], we 
see that \PV t * b V td \ = V tc , while |(T - P)V* b V ud \ = V uc . 

As we will see below, the theoretical input is a predic- 
tion of the quantity \T— P\. There are basically two com- 
peting models for calculating nonleptonic decays: QCD 
factorization (QCDf) [8(] and perturbative QCD (pQCD) 
0. An effective-theory approach to nonleptonic decays 
that incorporates the ideas of QCDf and pQCD is the 
soft collinear effective theory (SCET) [l0|. Here, we will 
only consider QCDf and pQCD. 

Many references have concentrated on the differences 
between QCDf and pQCD. However, here we focus on 
their similarities. For example, both are models of non- 
leptonic B decays where the long-distance and short- 
distance physics are factorized at the b scale. The short- 
distance physics is calculable as an expansion in the 
strong coupling constant, while the long-distance effects 
are encoded in the form factors, light-cone distributions, 



etc. In particular, both QCDf and pQCD have the same 
mathematical formula for T — P. As such, despite the 
different methods of calculation, the prediction of the 
central value and error on \T — P\ will be quite similar in 
QCDf and pQCD. We treat these models in turn in what 
follows. 

We begin with QCD factorization. QCDf finds that 



T 

P 
4! 



o a 4Biy + 03 + 2/?4 — -^PzEW 



2 a 4EW + 03 + 2/?4 — 2^3EW 



~ 04EW '(9) 

/W (io) 



where the normalization factor in the denominator is 
A a d = Ml„F^ d ^ K (0)f K G F /V2. The quantities con- 

tributingto T/A d and P/A° d are defined and computed 
in Refs. [iliH. 

In QCDf the various hadronic quantities are calcu- 
lated using a systematic expansion in 1/mf,. How- 
ever, a potential problem occurs because the higher-order 
power-suppressed hadronic effects contain some chirally- 
enhanced infrared (IR) divergences. In order to calculate 
these, one introduces an arbitrary IR cutoff. The key 
observation here is that the difference = T — P is free 
of these dangerous IR divergences. This difference has 
been calculated [13]. The latest value of |A d | = |T - P\, 
which is the only physical quantity, is: 



| A d | = (2.96 ± 0.97) x 1(T' GeV 



(11) 



Note that in the difference T—P, the effect of electroweak 
contributions has been neglected in this calculation, but 
their impact is expected to be very tiny. 

The fact that |A d | is known and is free of IR diver- 
gences was first used in Ref. [l2[ to predict the values of 
the various quantities in B s — > KK. Note that in a A/nib 
expansion, we can expect corrections of up to of O(10)% 
with respect to the leading term for a generous range 
of values of A and m&, although they are often found 
to be smaller than expected by dimensional arguments 
(for instance, see Ref. [H|]). However, this is not true if 
power-suppressed IR divergences arise. These can induce 
errors up to 100% (for example, see Ref. [IH). Although 
A d has a residual sensitivity to the finite and suppressed 
A/rrib corrections, being free from IR divergences makes 
it more robust. Other ways of trying to approach the 
problem of 1/rrib corrections can be found in Ref. 14 1. 

In pQCD, the calculation of various hadronic quanti- 
ties is different than in QCDf. In particular, all quantities 
are finite. Thus, the calculation of T and P in pQCD is 
different than in QCDf. However, the difference T — P 
is the same in both formulations. Indeed, in Ref. [l5| . 
it is noted that the mathematical expression for T — P 
is the same as in QCDf. In calculating this expression, 
the only difference is the renormalization scale, /i. QCDf 
takes (j, ~ nib [16], while pQCD uses fi 2 = A QOD mb- 
Although the difference T — P has not been computed 
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A* 


K 
Oil 






m s (2 GeV) 


mc/mb 


l(A d )| 


19.7% 


4.1%, 


0.9% 


15.8% 


5.7% 


53.7% 



TABLE I: Relative impact of each error on |A«j|, denned as 
of/^totoii with « = M> "f ! a 2, Fq^ k , fh s , m c /m b . 



explicitly in pQCD, the above suggests that its value is 
similar to that of QCDf. For this reason, in what follows 
we use the QCDf calculation and assume that |Ad| is as 
in Eq. |jnj. 

At present, the errors on the quantities in |Ad| are 
quite large so that the extraction of a (below) suffers 
significant theoretical uncertainty. There are four sources 
of error. They are (i) the renormalization scale /i, which 
enters a s (jj,), Ci(/i) [H and r^(fj) [13, (ii) the Gegen- 
bauer coefficients and which enter the light-cone 
distributions (their dependence on (i is irrelevant here 
given the large error taken for them) , (iii) the form factor 



^0 



11], and (iv) the quark masses m s (2 GeV) [l| and 
m c /mb. Note that the same sources of error also affect 
the pQCD calculation. In calculating the errors in |Ad|, 
these quantities are taken to lie in the following ranges 
0,113 : m b/ 2 < A* < 2m b, a i = 0.2±0.2, af = 0.1±0.3, 



^0 



*K 



0.34 ± 0.05, m s (2 GeV) = 103 ± 20 MeV, 
0.26 < m c /m b < 0.36. 

It should be noted that in pQCD, the form factors are 
calculable in terms of the light-cone distributions of the 
mesons 17| ■ In this approach it is assumed that the soft 
contributions to the form factors are suppressed. There is 
no rigorous estimate of the error in the form-factor pre- 
dictions from this assumption. In the QCDf approach 
the form factors are inputs which may be taken from 
QCD sum rule calculations [3]. The sum rule approach 
assumes that the soft contributions dominate the form- 
factor calculation. The inputs to this calculation then 
include, in addition to the light-cone distributions, ad- 
ditional parameters associated with the Borel transfor- 
mation. Here too there is no rigorous estimate of the 
error in the form-factor calculations from the assump- 
tions involved in this approach. The predictions of the 
form factors from the sum rule and the pQCD approach 
are similar and it is reasonable to guess that the error in 
the form-factor calculation in both approaches is similar. 
Rigorous information about the form factors may be ob- 
tained from lattice calculations in the future with further 
progress in this area. Alternately, the form factors may 
be measured in experiments where one can have access 
to B — ► K form factors in B — > Kll decays. One can 
also relate the B — » K form factors to B — > 7r form fac- 
tors, obtained from semileptonic decays, by using SU(3) 
symmetry if the SU(3)-breaking effects can be reliably 
estimated. 

Table 1 indicates the relative error associated with each 
of the inputs entering | A^|. The error related to /j, is quite 



large. However, this error will be reduced when higher- 
order calculations are done. Concerning the Gegenbauer 
coefficients, to be conservative we prefer to use the val- 
ues given in Ref. 



11] 



even if the errors are a bit too 
large. There have been several interesting developments 
on the lattice [lj| and with QCD sum rules [20| . leading 
to more precise numbers for these coefficients. However, 
the problem is always how to estimate the error associ- 
ated with the higher Gegenbauer modes not included in 
the truncated series. Hopefully, this difficulty will be re- 
solved in the future. The largest source of error is the 
ratio m c /mh. Clearly a greater control of m c will lead 
to a substantial reduction of this error. A better un- 
derstanding of nonperturbative effects may help in the 
extraction of m c from charmonium data [2l| . 

Now, both QCDf and pQCD claim that, at leading or- 
der in 1/rrib, the short- and long-distance physics in non- 
leptonic decays factorize. If nonfactorizable long-distance 
effects are important, the predicted value of |Ad| will 
change. Thus, any method of determining weak-phase in- 
formation that uses input from theoretical models based 
on factorization must also provide an independent test of 
such models which does not use that phase information. 
As we will see, the method described in this paper does 
just that. Thus, we propose a new way of obtaining a 
and a clean test of theoretical hadronic models based on 
factorization. 

We have shown that Vtc can be defined in terms of a 
and the experimental B ( j — » K°K° observables [Eq. flBJ)]. 
Vuc can be obtained similarly: 



V 



2 

uc 



z R -x 

cos 2a — 1 



(12) 



Now, using the relation \(T — P)V* b V u d\ = V U c, one can 
obtain the weak phase a. 

In order to extract a, one proceeds as follows. It is 
first necessary to relate the quantities X, Y, Z r , Z a to 
the measurements. It was noted in Ref. [22] that 



ns - /) 



Pc 



&irm 2 B 



\A(B^f)\\ 



BR(B — * /) = r(fl-/)5 



(13) 



where p c is the momentum of the final-state mesons in 
the rest frame of the B, and tub and To are the mass 
and lifetime of the B° meson, respectively. A similar 
expression holds for the averaged branching ratio, which 
includes both A and A. Thus, we have 



T0Pc 



(14) 



where B is the averaged branching ratio. For Y, we have 
Y 



X 



»(///' 



Y = kBA. 



dir 



(15) 
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Now recall that A m i X 
e -^A/A. Thus, 



-2Im(A)/(l + |A| 2 ), where A 



Z, = \A\ 2 lmX = -XA mix = -KBA mix . (16) 
Combining these expressions for X, Y and Z 1 , we get 



Z R = ±kB\ 1 



A 2 



A 2 



(17) 



Certain references adopt a different convention for A m i X 
which differs from that given above by a sign. However, 
as can be seen from the above expression for Z R , this sign 
is unimportant here. Eq. \12\ can now be solved for a: 



cos 2a 



- 2 sin a 



\v: b v ud \ 2 \A B 



1 ± 



~\ - A 2 -A 2 

dir mix 



(18) 



Note that Ref. [l2| also shows that there exists a relation 
(called a sum-rule) between | A<z|, the branching ratio and 
CP asymmetries of B® — > K°K°, and 7. The above 
equation is equivalent to this sum-rule. (Fo r a similar 
relation for the weak mixing angle (f> s see |23|). 

The various experimental quantities in B® — > K°K° 
have now been measured [24j |. They are 



B 

Adir 
Adir 

A 

- 1 111 / r 



0.95: 



0.19 



x 10" 



0.58l||±0.01 [Belle], 
-0.40 ± 0.41 ± 0.06 [Babar], 

-, 9O+0.73+0. 16 
-^"-O.SO-O.ll ■ 



(19) 



Combining these measurements (so far not very precise) 
with the value of \V* b V ud \ = (3.586tg;J?g) x !0~ 3 [H, 
and that for |A<j| given in Eq. (jlip . we can obtain sin 2 a. 
Unfortunately, at present, there is no constraint; all val- 
ues < sin 2 a < 1 are allowed. 

There is, however, more that can be learned. If one 
assumes that there is no new physics, then one can use 
other determinations of the weak phases. Along with the 
measurement of the B® — > K°K° branching ratio and the 
theoretical input of \A d \, knowledge of the weak phases 
allows one to make predictions for the B® — > K°K° CP 
asymmetries [Eq. (15} ]. At present, (3 and 7 are measured 
more accurately than a, so it is useful to write sin 2 a = 
sm 2 {/3 + 7). Taking (in degrees) (3 = 22M±° 7 6 2 2 and 7 = 
59.0IJ; 2 . (25)], the SM currently gives 0.95 < sin 2 a < 1. 
Now, the value of the coefficient (KB)/{\V* b V ud \ 2 \A d \ 2 ) 
is 98 ± 67. This requires A d i r and A m i X to be small, 
with the negative sign of the square root, in order for the 
right-hand side to reproduce the SM value of sin 2 a. In 



we obtain 0.02 < A% r - 



< 



particular, from Eq. (T 
0.125. This then is a prediction of this method. Put 
another way, this is a test of theoretical hadronic models 
based on factorization. If this constraint on A 2 dir + A 2 mix 
is not satisfied, it means that the expected SM value of 
sin 2 a is not reproduced, which in turn implies that the 
value of I Ad I [Eq. (fTTj) ] is wrong (or there is new physics). 



If the test fails, it means there are effects that QCDf 
and pQCD have not considered, i.e. nonfactorizable long- 
distance physics. One example of such an effect could be 
possible long-distance contributions from cc intermediate 
states to the c-quark piec e of the penguin diagram, P c 
("charming penguins" [lOLl^l. In this case, the calcula- 
tion of I Ad I in Eq. (flTj) is not correct, as it does not take 
these effects into account, and the test will fail. Though 
the impact of charming penguins on B decays is a con- 
troversial issue 27J, the method proposed in this paper 
for extracting a can be understood also as a test to find 
long-distance effects such as charming penguins. 

Note that if long-distance effects turn out to be impor- 
tant, the method proposed here for determining a can 
still be used if a reliable theoretical evaluation of T — P 
which includes the long-distance effects can be found. 

In summary, in this paper we have proposed a new 
method of obtaining the CP phase a. It involves mea- 
surements of the decay B® — > K°K°, along with the- 
oretical input. This input is largely the same whether 
it is taken from QCD factorization (QCDf) or perturba- 
tive QCD (pQCD). The theoretical error is very large at 
present. We have computed the percent errors due to the 
various contributions. We note that the theoretical error 
depends on quantities such as the renormalization scale, 
light-cone distributions, and form factors, which will all 
be better known in the future, so that this error will be 
correspondingly reduced. Unfortunately, the theoretical 
error also depends strongly on quark masses, particularly 
m c /m b . It appears difficult to increase the precision on 
this ratio, but if this can be done this method can even- 
tually provide a new way of extracting a, and may help 
to reveal the presence of new physics. 

If one assumes that new physics is not present, the 
method allows one to make predictions for the CP asym- 
metries in B® — » K°K° decays. This provides a test 
of the standard model, and permits one to distinguish 
between factorization models of nonleptonic B decays 
(such as QCDf and pQCD) and those which include long- 
distance effects. In this way it may be possible to ascer- 
tain which kind of physics is present in b decays. 
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